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ABSTRACT 

We present preliminary results from a deep (600 ks) Chandra observation of 
the hot interstellar medium of the nearby early- type galaxy Centaurus A (Cen A). 
We find a surface brightness discontinuity in the gas ~3.5 kpc from the nucleus 
spanning a 120° arc. The temperature of the gas is 0.60±0.05 and 0.68±0.10 keV, 
interior and exterior to the discontinuity, respectively. The elemental abundance 
is poorly constrained by the spectral fits, but if the abundance is constant across 
the discontinuity, there is a factor of 2.3±0.4 pressure jump across the discon- 
tinuity. This would imply that the gas is moving at 470±100 km s" 1 , or Mach 
1.0±0.2 (1.2±0.2) relative to the sound speed of the gas external (internal) to 
the discontinuity. Alternatively, pressure balance could be maintained if there is 
a large (factor of ~7) discontinuity in the elemental abundance. We suggest that 
the observed discontinuity is the result of non-hydrostatic motion of the gas core 
(i.e. sloshing) due to the recent merger. In this situation, both gas motions and 
abundance gradients are important in the visibility of the discontinuity. Cen A 
is in the late stages of merging with a small late-type galaxy, and a large discon- 
tinuity in density and abundance across a short distance demonstrates that the 
gas of the two galaxies remains poorly mixed even several hundred million years 
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after the merger. The pressure discontinuity may have had a profound influence 
on the temporal evolution of the kpc-scale jet. The jet could have decollimated 
crossing the discontinuity and thereby forming the northeast radio lobe. 

Subject headings: galaxies: individual (Centaurus A, NGC 5128) - X -rays: galax- 
ies - hydrodynamics - galaxies: jets 



1. Introduction 



Chandra has revealed the complex relationship between the lobes of radio galaxies and 
the hot gas of galaxies, groups, and clusters in a large number of systems. The heating of the 
gas induced by the inflation of the lobes, either subsonic or supersonic, plays an important 
role in the energy budget of the ISM of galaxies and the ICM of clusters. Lobe inflation also 
plays a key role in mixing t he low entropy gas in the cores with the high entropy gas in the 
halo (IRoediger et al. 1 120071 ). In addition, the relatively low inertia of the radio lobes means 
that they are easily influenced by subsonic or transonic motions of the ambien t medium once 
the early, highly supersonic expansion phase has ended ([Reynolds et al. Il200ll ). The temporal 
evolution of radio lobes could be dramatically altered in systems that have undergone recent 
mergers. Observationally, there are at least two examples in whi ch ongoing galaxy merg - 
ers have dramatically influenced the evolution of radio bubbles ( IHardcastle et al. I l2007bl ; 



Worrall et al. 120071). In Chan dra observations both of the Antennae, the nearest early- 



stage merger jBaldi et al. Ibood ). and of more distant merging galaxies (IBrassington et al. 



20061 ) we see large spatial variations of gas temperature, density and elemental abundance 
as the hot ISM of the merging galaxies begins to coalesce. The dissipative processes that 
smooth out these gradients are poorly understood, but transonic gas motions and temper- 
ature/pressure variations could play an important role in the time evolution of relatively 
weak nuclear outbursts. 

Centaurus A (Cen A, NGC 5128) is the nearest radio galaxy and the neares t large, un- 



obscured, early-type galaxy (D L =3.7 Mpc, M B =-21.2, see iDufour et al. I ( 119791 )). It is the 
dominant member of a poor group, and the prototypical example of several astrophysically 
interesting categories, including low-luminosity (FR I) radio galaxies, early-type galaxies, 
and late-stage mergers. These make it the ideal target for the study of the dynamics of 
the merger process and its effects on radio bubbles. Six 100 ks observations were made 



with Chandr a / AClS-l in AO-8 to study the X-ray jet (IKraft et al. 1 120021 : IHardcastle et al. 



2003 



2007b 



2006 



Worrall et al. 1 120081 ) . the interaction o f the radio lobes with the ambient 



gas (Kraft et al. 1120031). t h e X-ray binary popu lation (IKraft et al. I l200ll; IVoss fc Gi 



Jordan et al. 1 120071 : ISivakoff et al. 1 120081 ). and the hot gas (IKarovska et al. 



fanov 



2002 
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Kraft et al. 1120031 ) 



In this Letter, we present preliminary results from our deep AO-8 Chandra observation 
of the ISM of Cen A and its interaction with the northeast radio lobe. In particular, we 
report the discovery of a surface brightness discontinuity in the gas ~3.2 kpc from the 
nucleus that spans a 120° arc. We conclude that this discontinuity and associated gas 
motions may be responsible for the observed radio morphology of the northeast radio lobe. 
We adopt a distance of 3.7 Mpc to Cen A, which is the mean of distance measures from 
five indepen dent indicators (TRGB , Mira variables, SBF, PNLF, Cepheids, summarized in 



Section 6 of iFerrarese et al. I (120071 )). At this distance, 1"=17.9 pc and l'=1.076 kpc. All 



uncertainties are at 90% confidence for one parameter of interest unless otherwise stated, 
and all coordinates are J2000. All spectral fit s include absorption (iV#=8.41xl0 20 cm -2 ) by 



gas in our Galaxy (jDickey &: Lockman Ill990l ). 



2. Data Analysis 

Cen A has been observed ten times with Chandra/ ACIS. All the data were reprocessed 
to apply the most up to date gain and CTI correction, and the event files were filtered 
to remove events at node boundaries. Standard ASCA grade filtering (i.e. event grades 
0,2,3,4,6) was applied to the data. A lightcurve was created in the 5-10 keV band after 
removal of point sources to search for periods of background flaring. Intervals where the 
background rate was more than 3cr above the mean were removed, leaving ~623 ks and ~96 
ks of good time for the ACIS-I and ACIS-S observations, respectively. The observations 
were co-aligned by matching the X-ray positions of ~200 X-ray binaries and other point 
sources (about half h aving >100 counts in the 0.5-7.0 keV band within a single observation) 



( 1 Jordan et al. 1 120071 ). Full details of the data processing and alignment will be presented in 
a future publication. Given the narrow scope of this Letter, we use only the ACIS-I data 
(roughly 85% of the observation time) to avoid systematic uncertainties due to different 
backgrounds and spectral responses between the ACIS-S and ACIS-I detectors. 



3. Results 

A Gaussian smoothed, exposure corrected Chandra/ ACIS-I image of Cen A in the 0.5- 
1.0 keV band with 5 GHz radio contours (6" resolution) overlaid is shown in Figure [H One of 
the most striking features of Figured] is the surface brightness discontinuity which lies ~3.2' 
(3.54 kpc) to the east and north of the nucleus. Figure [2] contains an unsmoothed X-ray 
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image in the same band in the vicinity of the discontinuity (denoted by the blue arrows). 

We extracted a surface brightness profile, shown in Figure [31 across the discontinuity in 
a 50° pie slice region centered at a=13:25:34.445, £=-43:01: 16.84. The vertex of this region 
was located ~1.3' east of the nucleus, so that the radius of curvature of the discontinuity 
matches that of the annular regions, simplifying the deprojection of surface brightness to 
density. We fitted spectra to four radial regions (la, lb, 2a and 2b) of the galaxy gas. 
Background was determined from a distant region of the combined image free from emission 
from Cen A. Absorbed APEC models were fitted to each spectrum, with the absorption fixed 
at the Galactic value, and the abundance allowed to vary freely. The two spectra on the 
interior of the discontinuity (la and lb) contained an additional component, the contribution 
to the emission from the overlying gas. The best fit temperatures and 90% uncertainties 
are shown in Table [TJ In all cases, the elemental abundance was unconstrained. At these 
temperatures, the emission is line dominated so that the normalization can be traded off with 
the abundance. The only formal statistical constraint from each fit was that the elemental 
abundance, Z, is >0.1 at 90% confidence in the two regions interior to the discontinuity 
(regions la and lb). 



4. Interpretation 

We deprojected the surface brightness profile assuming a uniform gas density out to 
Rdiscont (134"), and a beta-model gas density profile with core radius (fixed) r =200" beyond 
Rdiscont centered at a point ~1.3' east of the nucleus (see Figure [2]). Our general conclusions 
are not sensitive to these model assumptions. The temperatures of the interior and exterior 
gas were taken to be 0.6 and 0.68 keV, respectively, and represent the average values of regions 
la+lb and 2a+2b. The elemental abundance was unconstrained in the spectral fits, so we 
initially assumed uniform abundance (Z=1.0) on both sides of the discontinuity. The best 
fit surface brightness profile using the density profile described above is shown overplotted 
on Figure [3J We find the hydrogen density interior and exterior to the discontinuity to 
be 3.60xl0" 3 and 1.38xl0" 3 cm , respectively. The density ratio is more than double the 
ratio of the temperature differences and implies that the interior gas is greatly overpressurized 
relative to the exterior gas. In this case, the interior gas should be expanding at roughly its 
internal sound speed and the pressure gradient would dissipate in roughly a sound crossing 
time. 



One possible interpretation of this phenomenon is that the surface brightness disconti- 
nuity and density jump are the result of 'sloshing', that is non-hydrostatic motio ns of the 
central dense regions of the ISM due to the merger (IMarkevitch &: Vikhlinin 1 120071 ). In this 
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scenario, the gas is moving to the northeast at transonic velocity and is oscillating, perhaps 
non-radially, in the gravitational potential of the galaxy. This model has been invoked to 
explain the presence of similar features observed in clusters of galaxies. Hydrodynamic sim- 
ulations of this phenomenon for galaxy clusters indicate that the non-hydrostatic oscillation 
of the core can be quite long lived, typically several Gyrs. The length scales are much smaller 
in Cen A than in galaxy clusters, but the sound speed of the gas is also much lower. The 
presence of the 'slo shing' gas is entirely consistent with the time since the merger (several 



hundred Myr - see llsrael I ( 119981 )). The ratio of the gas pressures interior and exterior to 



the discontinuity is 2.3±0.4. If this ratio represents the ratio of pressures at the stagnation 
point and free stream region, the gas interior to the discontinuity is moving to the northeast 
at approximately 470 km s^ 1 , or with Mach number 1.0±0.2 (L 2±0.2) relative to the gas 
exterior (interior) to the discontinuity (ILandau &: Lifshitz Ill950l ). 



In this analysis we assumed, however, that the elemental abundance is constant across 
the discontinuity. If there is a large, discontinuous jump in the elemental abundance, the 
gas on the two sides of the surface brightness discontinuity could be at or near pressure 
equilibrium. This would imply that the elemental abundance of the material interior to the 
discontinuity was several times larger than that exterior to the discontinuity, and that the 
ratio of the gas densities is the reciprocal of the temperature ratio. In this model, the surface 
brightness discontinuity is a relatively stable, long lasting contact discontinuity between two 
fluids. The cooling function, A, of an optically thin plasma scales nearly linearly with 
elemental abundance at the temperatures (<0.75 keV) and abundances relevant to the hot 
ISM of Cen A. The difference from a purely linear relationship is less than 5% for Z >0.1, 
and since uh ~ A -0 ' 5 , these differences are negligible. The elemental abundances of the gas 
interior and exterior to the discontinuity would therefore have to differ by a factor of ~7 for 
the two fluids to be in pressure equilibrium. Thus, if the exterior gas is relatively metal poor 
(Z ext =0. 1-0.3), the gas interior to the discontinuity would have to have Solar or super-Solar 
abundance (Z int =0. 7-2.1) for the pressure to be continuous. 

Such a low abundance of the gas exterior to the edge is plausible. Recent HI observations 
of a sample of relatively isolated elliptical an d lenticular galaxies d emonstrates that there is 
considerable cool gas residing in their halos (lOosterloo et al. 1 120071 ). If such primordial cold 
gas is continuously falling into Cen A, the elemental abundance of the outer parts of the 
hot ISM could be quite low. Arecibo measuremen ts of HI around early- type galax ies in and 
around the Virgo cluster have found little cold gas (jdi Serego Alighieri et al. 1120071 ) . In dense 
environments, HI is li kely to be st ripped, but this is not expected in a poor environment 
like the Cen A group (llsrael Ill998l ). A recent X-ray stu dy of the hot ISM of the gas poor, 



isolated early-type galaxy NGC 4697 found Z ~0.25Z Q (llrwin et al. 1 120071 ). consistent with 



this picture. There have been several measurements of the elemental abundance of the stellar 
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halo population of Cen A. lHarris fc Harris I (120001 ) reported that the average abundance 
of the stars in the halo of Cen A (~20-30 kp c from the nucle u s) is Z ~O.37Z . After 
subtraction of the large-scale halo contribution, lHarris fc Harris I (120021 ) found that the the 
metallicity of stars closer to the nucleus (~8 kpc) was somewhat more metal rich, with 
Z ~0.65Z Q . Suzaku measurements of emission lines and elemental absorption edges of the 
circumnu clear gas in Cen A sug gest that the elemental abundance is slightly super-Solar 
(Z=1.3) (IMarkowitz et al. 1 120071 ). However, Suzaku does not have the angular resolution 
to resolve any of the substructure in the gas, the variable absorption, or the X-ray binary 
population. Analysis of Chandra spectra of the gas in the inner 2 kpc of Cen A demonstrates 
that the gas temperature of the central region var ies by more than a factor of two, and there 
is a spatial variation in the elemental abundance (IBrassington II2006I ). The temperature and 
elemental abundance of the central region based on the VLP observations will be presented 
in a future publication (Kraft et al. 2008, in preparation). 

In the absence of a precise measurement of the elemental abundances across the dis- 
continuity, we conclude that both phenomena (i.e. gas motions and elemental abundance 
gradients) may be contributing to the visibility of the surface brightness discontinuity. The 
abundance profile in many early-type galaxies is sharply peaked at the center, and if the 
central core is rapidly displaced off-center, a sharp surface brightness discontinuity could 
result even at zero velocity. The inferred velocity of ~470 km s _1 should therefore be 
regarded as an upper limit. The effect of projection is another source of systematic uncer- 
tainty in this analysis. Three d imensional simulations of sloshing in galaxy cluster mergers 
( IMarkevitch fc Vikhlinin 1 120071 ) show complex gas temperature and density variations that 
would be difficult, if not impossible, to resolve via deprojection of surface brightness profiles. 

We consider the possibility that this surface brightness discontinuity is a shock due to 
the inflation of radio lobe to be unlikely for several reasons. First, it is not at all clear what 
is driving the shock. The pressure is larger on the interior, so the shock must be driven from 
the inside out. The northeast radio lobe clearly lies outside the discontinuity, so the inflation 
of the lobe cannot be responsible. There is nothing that appears to be driving the gas at 
supersonic velocity more or less radially from the nucleus. Second, the surface brightness 
and temperature profiles are wrong for a shock. We should observe a sharp jump in the 
surface brightness and gas temperature just behind the shock. This temperature jump, even 
for a high Mach number shock, would easily be observable by Chandra in Cen A if present. 
Third, there is no corresponding feature on the opposite side of Cen A. It is hard to imagine 
how the gas could be shock heated on only one side of the galaxy. 
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5. Implications for Evolution of the Jet and Radio Lobes 

Did the discontinuity in the gas affect the temporal evolution of the jet and northeast 
radio lobe? The X-ray morphology of the jet shows a nearly constant open ing; angle to a dis- 



tance o f 139" from the nucleus, where the X-ray jet narrows (see Figure 1 of lHardcastle et al 



( I2007bl )). It is at about this point that the jet opens into the northeast lobe in the radio 
band. One possibility is that the point at which the X-ray jet narrows is the position at 
which the jet crosses the X-ray surface brightness discontinuity. This point lies interior to 
the discontinuity in projection since the jet does not lie directly in the plane of the sky. The 
sudden change in ambient gas density and pressure across the discontinuity could explain 



the transition from jet to lobe. iNorman et al. I (119881 ) modeled this jet to lobe transition in 



Cen A as the result of the jet encountering a sudden increase in the external pressure and 
the jet transitioning from supersonic to subsonic relative to its internal sound speed. They 
hypothesized that gas was infalling on Cen A and driving a shock in toward the center the 
galaxy. Our data suggest that it is, in fact, the opposite scenario that is occuring. The 
motion of the gas interior to the discontinuity is creating the pressure jump (decrease) that 
de-collimates the jet. 

The decollimation of the jet across the discontinuity may account for some of the ob- 
served features of the X-ray jet. There must clearly be a change in the efficiency of particle 
acceleration at the point where the jet cross the disc ontinuity. The X-ray em ission from the 



jet narrows and fades away as the jet enters the lobe (lHardcastle et al. Il2007bl ). The spectral 
index of the jet emission beyond the discontinuity is considerably steeper than either the 
compact knots or diffuse emission of the rest of the jet. There is no obvious structure in the 
radio emission, however. The radio jet both widens and brightens as it smoothly enters the 
northeast radio lobe. Rapid expansion of the jet may lead to a reconfinement shock, causing 
the jet to decelerate, and will certainly cause the magnetic field strength in the jet to fall: 
both of these could reduce the efficiency of acceleration of ultra-relativistic particles. 

Can the change in physical conditions at the discontinuity be related to the bending 
of the lobe toward the northwest? Fundamentally, the bending of the X-ray jet into the 
northeast lobe is likely the result of velocity and pressure gradients in the external gas, 
although the low surface brightness of the gas and uncertainties due to projection prevent us 
from making a definitive statement. If the core is 'sloshing' there will be significant subsonic 
or transonic motions of the gas exterior to the discontinuity. This implies pressure variations 
of tens of percent to perhaps as large as a factor of 2. However, sharp X-ray surface brightness 
discontinuities have been observed along the north and northwest boundaries of the northeast 



lobe suggesting that it is expanding into the ISM supersonically (IKraft et al. 1120071 ). We also 



know that the southwest radio lobe is expanding into the ISM highly supersonically. The 
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similarity in size of the northeast and southwest lobes demonstrates that they must be 
inflating at roughly equal rates, although there are obvious bends in the jet in projection, 
so there may be some difference in the sizes of the lobes. A nearly transonic wind striking a 
radio lobe could cause it to inflate asymmetrically if the lobe is overpressurized by at most 
a factor of a few relative to the ambient medium. 

Finally we note that the position of the X-ray surface brightness discontin uity perfectly 



overl a ps the position of one of the optica l shells seen in deep exposures of Cen A (IMalin et al. 



19831 ; iGopal-Krishna fc Saripalli Ill984l ). These optical shells (in Cen A and other early- 
type galaxies) are believed to be the result of phase wrapping of the dynamically cold disk 
of the merging spiral galaxy in the gravitational potential of the more massive early-type 
galaxy (jQuinn 1 11984 ) . The motions of the stars in these shells are thought to be a small 
perturbation in the larger scale potential. If the association of the X-ray surface brightness 
discontinuity with the stellar shell is not coincidental, the gas and some of the associated 
stars must therefore be moving synchronously. Interestingly, both H I and CO have been 



observed to be associated with the optical shells of Cen A as well (jSchiminiovich et al. 



1994 ; ICharmandaris. Combes, fc van der Hulst 1 120001 ) . The appearance of the optical shells 
may not be the result of phase-wrapping but could be indicative of other dynamic stellar 
motions. 
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Fig. 1. — Gaussian smoothed, exposure corrected Chandra/ ACIS-I image of Cen A in the 
0.5-1.0 keV band with 5 GHz radio contours (6" resolution) overlaid. The color stretch has 
been selected to enhance the appearance of the X-ray surface brightness discontinuity to the 
northeast of the nucleus. 
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Fig. 2. — Raw Chandra/ ACIS-I image of Cen A in the 0.5-1.0 keV band. The light blue 
arrows denote the approximate position of the surface brightness discontinuity in the gas. 
The light blue cross and yellow lines denote the wedge in which the surface brightness profile 
in Figure [3] was created. The approximate position of the transition region of the jet (where 
the width of the X-ray jet narrows) is also labeled. 
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Fig. 3. — Surface brightness profile of the hot ISM in a 50° wedge (shown in Figure [2]) in 
the 0.5-1.0 keV band. The continuous curve is the best fit surface brightness profile model 
consisting of a uniform density sphere out to radius 134" and a beta model profile at larger 
radii. 
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Region 


Temperature (keV) 


la 


0.61±0.05 


lb 


0.59±0.04 


2a 


0.73±0.10 


2b 


0.63±0.09 



Table 1: Best fit temperatures and uncertainties (90% confidence) for two regions interior 
(la and lb) and two regions exterior (2a and 2b) to the surface brightness discontinuity. 



